Is the ultra-fast transformation of bismuth non-thermal? 

E. G. Gamaly, A. V. Rode 

Laser Physics Centre, Research School of Physics and Engineering 
Australian National University, Canberra, ACT0200, Australia 

Abstract 

Transient state of femtosecond laser excited bismuth has been studied by various groups with 
time-resolved optical, x-ray, and electron probes at the deposited energy density from below 
through up to several times the equilibrium enthalpy of melting. However, the interpretations of 
the experimental results are controversial: the optical probes reveal the absence of transition to 
the melting phase while the authors of x-ray and electron diffraction experiments claim the 
observation of ultrafast non-thermal melting. The presented analysis, based on temperature 
dependence of bismuth optical properties, unequivocally shows a purely thermal nature of all the 
observed fs-laser induced transformations in bismuth. 

I. INTRODUCTION 

The advent of ultra-fast optical, x-ray, and electronic diagnostic techniques have provided an 
insight into such most common phase transformation as crystal melting at atomic space and time 
scales. The transient state of the femtosecond laser excited Bi crystal has been studied with time- 
resolved fast optical [1-10], x-ray [11-13] and electron [14,15] probes at the deposited energy 
density below and up to eight times above the equilibrium enthalpy of melting. However, there 
is no consensus in interpretation of the experiments with the same crystal in nearly identical 
excitation conditions. The fast fall-off of x-ray probe intensity reflected from bismuth crystal 
[12] as well as the intensity decrease for the diffracted electron beam [14] were interpreted as the 
evidence of non-thermal lattice melting induced solely by the electronic excitation, while the 
dual optical probe of bismuth reflectivity in the same conditions clearly demonstrate the absence 
of transition to melting state and presents an evidence of purely thermal character of the ultra-fast 
material transformation [7,8]. 

A recent publication of Sciaini et al. [14] is another striking justification of the existing 
controversy. The authors have examined the transient states of the 30-nm freestanding bismuth 
films irradiated by 200-fs laser pulses through femtosecond electron diffraction. It is remarkable 
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that the sharp decrease in the diffracted intensity was observed even at the absorbed energy 
density half of the equilibrium melting enthalpy when the disordering due to melting is 
energetically impossible. Obviously, the reason for the observed intensity decrease calls for a 
different explanation. 

The measured temperature dependence of optical properties of bismuth in equilibrium 
conditions in the temperature range from room temperature up to 773 K (-21 5K over the melting 
point) [7,10,17-19] allows one obtaining the high temperature properties of Bi directly from the 
experiments. The material properties recovered from the optical data, such as the electron- 
phonon momentum and the energy transfer rates, the heat diffusion coefficient, the number of 
electrons in the conduction band and others, form a strong basis for the quantitative interpretation 
of the experimental results with the ultra-fast excitation of bismuth. The time for setting up the 
statistical distributions and the electron-lattice thermalisation time also directly follow from the 
optical data. The diffusion coefficient recovered from the optical data is in close proximity to 
that obtained from x-ray reflectivity in fs-excited bismuth experiments by Johnson et. al. [13]. 

The structure of the paper is as the following: (i) - we retrieve the temperature dependencies 
of the electron-phonon momentum transfer rate and the plasma frequency (the number of 
electrons in the conduction band) obtained from the dielectric function dependency on 
temperature in the equilibrium conditions using the Drude-like form for the dielectric function of 
bismuth; (ii) - we derive theoretically the temperature-dependent electron-phonon momentum 
exchange rate and its link to the energy exchange rate. We fit the momentum rate to the 
experimental data and then establish the temperature-dependent energy exchange rate 
corresponding to the experiments; (iii) - we estimate the absorbed laser energy and the lattice 
temperature in Bismuth layer excited by laser with fluences (0.5 - 23) mJ/cm ; (iv) - we compare 
the electron-phonon energy transfer time to the characteristic time of decrease in the electron and 
x-ray diffracted intensity to the quasi-stationary level as the function of the absorbed laser 
fluence; (v) - we identify the processes, which may induce non-homogeneous heating of the 
skin-layer, which, in turn, are the reasons for the breaking of the Bragg conditions and therefore 
for the observed diffracted intensity decrease. Finally, we draw the conclusion on the nature of 
transformation in bismuth excited by the fs-laser pulse. The analysis presented in this paper 
unequivocally shows a purely thermal nature of all observed transformations in bismuth excited 
by ultra-fast laser pulses. 
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II. TEMPERATURE-DEPENDENT BISMUTH PROPERTIES FROM OPTICAL 
EXPERIMENTS IN EQUILIBRIUM 

A. Electron-phonon momentum exchange rate 

The measurements of the reflectivity and the dielectric function by the ellipsometry technique, 
which gives simultaneously the real and imaginary parts of dielectric function allow direct 
recovering the electron density and the electron-phonon momentum exchange rate. Indeed, the 
reflectivity, R, is directly related to the real e r and imaginary e t parts of the dielectric function, 
s= s r + i8 i through the Fresnel formula: 
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where ^ = ~\i r +s] . It was found long ago [17-20] that the dielectric function for solid Bi at 

room temperature and that for the liquid Bi obeys to the Drude-like form. The dielectric function 
in the Drude form is directly linked to the plasma frequency co p and the electron-phonon 
momentum exchange rate, v e . ph : 
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Here co is the laser field frequency, co = 2.356x10 s" for 800 nm, and co =Ane — e -, e is the 

m e 

electron charge and m* e is the effective electron mass. Thus the electron-phonon momentum 
exchange rate and plasma frequency, measured in units of laser frequency, are directly connected 
to the real and imaginary parts of the dielectric function: 
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If the electron mass is known than the number density of the conductivity electrons is directly 
retrieved from the plasma frequency. Following this procedure Garl [7] recovered the following 
Bi properties from the ellipsometry measurements at the room temperature T RT = 294K for laser 

2 

light of 800 nm: % = 31; R = 0.74; v e . ph = 2.1xl0 15 s" 1 ; n e = 5.34xl0 22 cm" 3 ; Fermi energy e F = 
co 

5.17 eV (vf= 1.35x10 cm/s) under assumption that electron has a free electron mass. Comins 
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[17] found that Bi at 773 K has the following parameters: The reflectivity is R = 0.67; (co pe /co) = 
81.58 that corresponds to the number density of free carriers n e = 1.42x10 cm" . Therefore all 5 
valence electrons are transferred into conduction band. Comins suggested that electron has a free 
mass value [17], then the Fermi energy ef = 9.92 eV, vf= 1 .87 x 10 cm/s in molten Bi. The 
electron-phonon momentum exchange rate reads v e . p h = 5.67><10 15 s" 1 . Therefore, the available 
data from 293 K to 773 K changing seemingly smoothly. However, there is a phase transition 
point from solid to liquid state, the melting point, at 544.7 K. It is known that Bi, along with ice, 
Sb, Ga, KI, and InSb, belongs to so-called open crystalline structure [21] so its molten phase 
density is higher than the solid state. Little is known about the behaviour of other physical 
parameters of Bi near the phase transition point from the literature. Moreover, there is no 
consensus, to the best of our knowledge, if the melting of Bi is the first order or the second order 
phase transition. As it follows from the text below, physical parameters recovered from the 
optical experiments such as number density of electrons in the conduction band and the Fermi 
energy value, also apparently smoothly changing with the temperature. We have no data near 
and above the melting point, so we decided to interpolate the available data by smooth 
continuous dependence on the temperature as a first approximation until some new data become 
available. Some support to this action comes also from the theory presented further down in the 
paper, where the numerical coefficients in the linear temperature dependence for the electron- 
phonon momentum exchange rate coincide with the experimental values and are almost identical 
for both solid and liquid states. Therefore, linear interpolation of the data extracted from all 
available experiments [7,17,18], and keeping in mind the theory (see in Section III.A below) 
shows that v e . p h grows up in direct proportion to temperature (Fig.l): 

^ A =2.1xl0 15 ^[5- 1 ]. (4) 

The linear dependence of the momentum exchange rate on temperature fits the optical 
measurements for bismuth in equilibrium with sufficient accuracy; the proportionality to the 
temperature holds for Bi well before and long after the melting point [7,17-19]. 
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Fig.l. Temperature dependence of the momentum exchange frequency (a), and of the electron 
density (b), both in equilibrium conditions; solid line corresponds to linear dependence 
Eq.(3,4). The circles are the results from ellipsometry measurements Ref.[7], triangle - 
Ref.[17], square -Ref.[ 18]. 



B. Number of electrons in the conduction band 

The electrons number density in the conduction band is retrieved from the plasma frequency 
under assumption that the electron mass is equal to a free electron mass: m = m e . The result at 

22 3 

room temperature gives n e = 5.34x10 cm" (1.89 out of total 5 valence electrons are in the 
conduction band) and e F = 5.17 eV (Fig. lb). At 773 K in liquid bismuth all five electrons are in 

23 3 

the conduction band: n e = 1.42x10 cm" , e F = 9.92 eV, and the dielectric function obeys the 
Drude-like form [10,17,18]. Thus the experiments had shown that from 40 to 100% of the 
valence electrons are transferred to the conduction band at the temperature increase from the 
room temperature to the melting point. These data are in a sharp contrast to the statements in 
[11,12,14] that only from 5 to 15% of valence electrons are in the conduction band under the 
strong laser excitation to a much higher temperature. 

C. Electronic heat conduction and characteristic cooling time 

Thermal diffusivity (coefficient of thermal diffusion) relates to both the Fermi energy and the 
electron-phonon momentum exchange rate: D -v 2 F l3v e _ ph . The diffusivity values recovered 
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from the temperature dependencies in equilibrium [10,17-19], are D = (2.0 - 2.89) cm /s for the 
temperature range 293 K - 793 K. These results are in good agreement with the recent non- 
equilibrium measurements from the x-ray reflectivity data of fs-laser excited bismuth giving D = 

2 2 

2.3 cm /s [13]. Note the drastic difference of this value from 0.067 cm /s given for equilibrium 

conditions in the reference book [22]. Accordingly, the electron mean free path, l mf = Vp = 

v , 

e- ph 

0.67 nm « l s , is much less than the skin layer depth and the film thickness, assuring the 
legitimacy of the diffusion approximation for the electron heat transfer. 

The time for the temperature smoothing across a 30 nm thick film [14] taken in accord to the 
diffusivity value is, respectively, t smooth =l 2 s I D =3.9 ps, which is a reason for preserving a high 
thermal gradient across the film thickness in the first ps after the laser excitation. 



III. LINK BETWEEN TEMPERATURE-DEPENDENT ELECTRON-PHONON 
MOMENTUM AND ENERGY EXCHANGE RATES: THEORY 

A. Electron-phonon momentum exchange rate 

The electrons are travelling in a periodic self-consistent potential W{r) of an ideal crystalline 
lattice without losses. The electron-phonon scattering is the main source of changes in the 
electron momentum and energy. The potential slightly perturbed (W » SW (r)) by the electron- 
phonon scattering, this could be presented as an expansion of W{r) into series of small 
perturbations [23]: 

W(r) = W + 5W{r) + 5W"(r) + ... = H Q + H' e _ ph {r) + H" e _ ph (r) + ... (5) 

The first correction to the potential corresponds to the emission or absorption of a phonon with 
quasi-momentum to be conserved. The perturbation of the potential expresses through the 
amplitude of atomic vibrations (phonons amplitude) A, as the following: 



or 



A, (6) 



The summation over the all phonon states is assumed. With the help of quantum-mechanical 
perturbation theory one can find the probability of electron transition from a state marked by 
indices, (p,l), to a state (p ',/'), simultaneously absorbing or emitting a phonon ha> kj in the kj 
mode [23]: 
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Here the delta-function represents the energy conservation. On the classic (as oppose to 
quantum) language the above probability is the electron-phonon momentum exchange rate. 

dW J A 

Let us estimate the perturbation Hamiltonian as H' = A = — — ; the value of the potential 

dr d 

is of the order of the first ionisation potential Jo, W=Jq (Jo = 7.3 eV for bismuth), while the space 
scale of the potential gradient is of the order of magnitude of the inter-atomic distance d. The 
electron energy is significantly larger than the phonon energy hco ph and it is practically 

unchanged in a collision, s F » ho ph . Then, Eq.(7) converts into the general expression for the 

electron-phonon momentum exchange rate v m e °™ h : 

At the temperature below the Debye temperature, T D , the electron-phonon momentum 
exchange rate (and the phonon's amplitude) does not depend on the lattice temperature, (T D = 
1 19 K for bismuth) k B T < hc s k « hco D , here c s is the speed of sound, k is the wave vector, and cod 
is the Debye frequency [23-25]. The phonon amplitude is A 2 «2h/ Mco D . Then, taking into 
account that Mco 2 D d 2 « s b , one obtains the collision frequency to be proportional to the Debye 

j 2 



frequency, v™™ h = <\n — —£0 D , at low lattice temperature. 
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At the higher temperature T L > T D one should average the squared phonon amplitude (energy) 
in Eq.(8) over the Maxwell-Boltzmann distribution. For ID oscillator this average expresses as 
the following [26,27]: 



} A 2 exp(-M« 2 A 2 /2k B T L yA 
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For harmonic oscillator in three dimensions the average squared phonon's amplitude reads: 

(9) 

1 Mco ph 
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The phonon frequency conventionally expressed through the second derivative of the inter- 
atomic potential in equilibrium: 



Mco 2 ph 



\dr 2 j 
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Sb = 2.16 eV for bismuth. Now Eq.(8) with the help of Eqs.(9,10) reads: 

T 1 k T 
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For Bi: Jo = 7.3 eV; sf= 5.17 eV, therefore the material-dependent coefficient in Eq.(l 1) 2n 

e F e b 

= 29.98, while the numerical coefficient a = (1-3) can be used for fitting to the experimental data. 
For the room temperature 294 K one gets for solid Bi from Eq.(ll) v™p h (T RT ) = 

(1.15-3.46) x 10 15 s" 1 . Thus, theory gives this rate proportional to the temperature in agreement 
with the optical ellipsometry measurements [7,17-20]. It is interesting that the experimental 
value of v e -ph = 2.1xl0 15 s" 1 (v e _ ph /co = 0.89) at room temperature could be attributed to the 
electron interaction with a two-dimensional harmonic oscillator at a = 2. 



B. Electron-phonon energy exchange 

The electron-phonon energy exchange includes the processes of emitting and absorbing 
phonons. Therefore the second term in the expansion of the interaction Hamiltonian should be 
taken into account: 



Hl ph =5W\r)* X - 



K dr 2 j 
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The electron-phonon energy exchange rate now expresses as the following: 

F 

The electron-phonon energy exchange rate (the inverse of the electron - lattice temperature 
equilibration time) is linked to the electron-phonon momentum exchange rate Eq.(8) as the 
following: 

A 2 

y en mom^_ ^ 

e-ph - e-ph V > 
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At the temperature below T D it reads: 



e— ph 
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(15) 



That is in a qualitative agreement to [28]. The Eq.(14) provides a basis for obtaining the energy 
exchange rate from the electron-phonon momentum exchange rate directly extracted from the 
optical experiments. One shall take the ratio of the squared phonon's amplitude to the inter- 
atomic distance in Eq.(14) corresponding to the two-dimensional oscillator as suggested by the 

experimental data (A 2 )/ d 2 — 2k B T L /s h . Then Eq.(14) takes the form: 
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Now, taking the link between the temperature dependent energy exchange and the collision 
frequency and using the experimental temperature dependence Eq.(4) one obtains the electron- 
phonon energy exchange rate dependence on temperature as the following: 
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Fig.2. Energy exchange rate: solid line corresponds to Eq.(17). Circles are the results from 
ellipsometry measurements Ref.[7], triangle - Ref.[17], square - Ref.[18] 
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The results of calculations of the energy exchange rate (with correction taken from 
experiments) are compared in Fig.2 to the available experimental data. Now it is also possible 
calculating the dielectric function and reflectivity changes with temperature, which are presented 
in Fig.3. The dielectric function at room temperature was measured to be s r = -16.25, = 15.4 
[7,10] at 800 nm; the values are very close to the literature data [20,22]. The electron-phonon 
collision rate and plasma frequency at room temperature were found with the help of Eqs.(2,3) as 
v e°ph = 0.893a* and co 2 pl = 3 1 .0ct> 2 . The temperature dependencies were calculated with the help 
of Eqs. (11,17). The theory predicts well the experimental results of dielectric function 
measurements for liquid bismuth at 610 K [18] and 773 K [17]. 
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Fig.3. Temperature dependencies of real s r and imaginary parts of the dielectric function (left) 
and reflectivity (right) of bismuth at 775 nm in equilibrium conditions. The solid lines 
were calculated using Eqs. (1,2, 17); the circles are results of ellipsometry measurements 
from Ref.[7], triangles - Ref.[17], squares - Ref.[18], arrows indicate the melting point 
T m = 544.7K. 



A. Absorbed laser energy, temperature and heat conduction 

The absorbed laser energy density E a b s in a skin layer excited by a fs-laser pulse is defined by 
the following relation [29]: 





Temperature (K) 



IV. ULTRAFAST LASER INTERACTION WITH BISMUTH 




(18) 
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where I(t) is laser intensity during the pulse of duration t p . A laser beam with the wavelength 775 
nm excites Bismuth in the skin layer l s = 28 nm with the absorption coefficient A = 0.26 [9,10]. 

The laser-excited Bismuth was studied in the fluence range, F(t p y (0.5 - 23) mJ/cm 2 [9-12] that 

"5 "5 

corresponds to the absorbed energy density range from 0.13 kJ/cm to 4.3 kJ/cm , i.e. from three 
times below through to eight times higher the equilibrium enthalpy of melting of 0.53 kJ/cm 
(that corresponds to the incident fluence of 2.85 mJ/cm ). The statistical distribution in a sub- 
system of identical particles is much faster than the laser pulse, providing a basis for the two- 
temperature description of the laser-excited matter [30,31,28]. The electron-electron collisions 
are responsible for establishing distributions in electron sub-system (electron temperature) [26], 
while the phonon-phonon collisions lead to the establishing the lattice temperature [23]. The 
absorbed energy density (intensity) and therefore the temperature exponentially decay inside the 
skin layer as the follows: 

T(x,t)=T maK (ty W {-2x/l s } (19) 

The temperature has a maximum at the vacuum-sample interface. The electron temperature 
decreases when electrons transfer their energy to the lattice via electron-phonon collisions and 
due to heat conduction later in time. As it follows from Fig.2 the energy exchange time for the 
studied range of laser fluences is shorter than the pulse duration. Therefore the electron-lattice 
equilibration occurs early in the pulse time and temperature (T e =T L ) reaches maximum at the end 
of the pulse. The maximum lattice temperature after the electron-lattice energy equilibration 
explicitly links to the absorbed laser energy due to the energy conservation: 



T L = ^ Eahs „ v ( 2 °) 



(C L n a + C e n e ) 



where Cl, C e are the lattice and electron heat capacity correspondingly, and n a , n e are the atomic 

n 2 k T 

and electron density in solid Bi [9,10]. We assume here that C L = 3k B and C e = — -, k B = 

1.38x10" J/K is the Boltzmann constant. The maximum temperature dependencies on the laser 
fluence are presented in Fig.4. 



The electron temperature at the vacuum-sample interface of the film is much lower than the 
Fermi energy s F = 5.17 eV even for the maximum fluence of 23 mJ/cm . Therefore the Fermi 
electrons with velocity v F = 1.35x10 cm/s are responsible for the heat transport across the film. 
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Thus, the temperature at the outer film boundary decreases due to the linear electronic heat 
conduction in accord to: 



T m ^{x = 0;t)~ T max (x = 0;t co J 



1/2 



(21) 



The cooling time for skin layer with the thickness l s is introduced above in the conventional 
form: 



cool s 



(22) 



here D is the coefficient of thermal diffusion, which is linked directly to the electron velocity, v F 
(the Fermi velocity), and electron-phonon momentum exchange rate, v e -ph, by the relation: 



D = v 2 F /3v , 



(23) 
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Fig. 4. Maximum lattice temperature after the electron-lattice energy equilibration vs incident 
laser fluence. The inset shows typical electron and lattice temperature behaviour during 
the laser pulse in the case when the energy equilibration time is shorter than the fs-laser 
pulse. 



B. Energy exchange time vs diffraction intensity fall-off in fs-laser excited bismuth 

The dependence of the energy exchange time f"_ ph = (v e Jl P hJ f° r me broad range of fluences 
in the laser-excited Bi is presented in Fig. 5 and compared to the measured characteristic time for 
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the diffraction intensity decrease in the experiments [11,12,14]. Analysis of the data from Fig. 5 
allows drawing at least two conclusions. First, it is remarkable that the intensity of the diffracted 
beam decreases significantly even at the fluences when the absorbed energy density is lower than 
the equilibrium enthalpy of melting. Therefore the melting is not the reason for the violation of 
the Bragg conditions and the intensity decrease. The laser-beam induced non-homogeneities in 
the excited sample might be responsible for the observed decrease. We discuss those later in the 
paper. The second conclusion from the analysis of Fig. 5 is that the times for the diffracted 
intensity decrease are always much longer than the electron-phonon energy exchange time, 
which represents the lattice heating time. Hence it follows, that the observed intensity decay is 
thermal in nature. 
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Fig.5. Electron-lattice energy exchange time t e e " ph = (v e e 'l P hJ according to Eq.(17) (solid line) 

and characteristic decrease time of the diffraction intensities: circles are from the electron 
diffraction experiments Ref.[14]; diamond (Ref.[ll]), triangle (Ref.[12]), and square 
(Ref.[13]) are from x-ray diffraction experiments. 

At the low fluence range the absorbed laser energy is up to five times lower than the enthalpy 
of melting. Therefore melting could not be the reason for the fall off of the diffraction intensity. 
The energy transfer from electrons to the lattice and the diffraction intensity fall-off occurs 
simultaneously at low intensity, while at superheating above the enthalpy of melting the 
thermalization process is much faster. Thus all processes at the absorbed energy density at and 
above the enthalpy of melting occur in conditions when electron and lattice temperature are 
equilibrated and the main part of the Maxwell-Boltzmann distribution (temperature) is already 
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established. Therefore all the transformations occurred are of the thermal nature. We discuss 
below the spatial anisotropy induced by the laser beam in the skin layer that might be a reason 
for the violation of the Bragg 's conditions. 



C. Laser-induced spatial inhomogeneity across the laser-excited film 

Let us discuss the possible sources for the laser-induced spatial inhomogeneities in the laser- 
excited sample ignored in the previous analyses. The spatial distribution of the absorbed laser 
energy across the skin layer might be the primary source of the observed fast drop of the 
diffracted beam intensity. Indeed, the absorbed laser energy decreases exponentially in the l s = 
28 nm bismuth skin layer: is afa (x)=.E' 0&smax exp^-2x/ / s ) - see Fig. 6. Therefore, both the absorbed 

energy density and temperature have maximum at the vacuum-sample interface irradiated by the 
laser pulse. The surface atoms start to expand after the energy transfer from the electrons with 
the speed of sound in Bi of 1.79xl0 5 cm/s. Thus in a picosecond time the outer interface moves 
on 17.9 Angstroms, more than 4 inter-atomic distances. Note that the gradient force 
F = —V(T e + T L )=(T e + T L yi skjn acts on atoms inside the skin layer slightly moving them in 

direction of the laser beam (into the sample) - Fig. 6. Definitely a coherent displacement affects 
the probe beam diffraction. There is no obvious way of distinguishing the decrease in the 
diffraction beam intensity caused by the inhomogeneity from that in the thermally disordered 
media. 

The laser interaction with a freestanding layer in [14] imposes further implications for 
interpretation of the electron diffraction pattern. The pressure at the front surface of the film, 
which is equal to the absorbed energy density, is ~35 kBar at 23 mJ/cm , while at the rear surface 
it is ~4 kBar. The sound velocity in solid Bi is 1.79 x 10 5 cm/s; thus the front surface of the free 
standing film shifts towards the laser beam during the 200 fs pulse by 3.6 A while the rear 
surface expands in the opposite direction. This creates additional source for the intensity 
decrease in the diffracted electron beam irrelevant to melting. Note that the temperature 
smoothing time due to heat conduction across the skin layer (see the end of Section II. C) is 3.9 
ps, which is much longer than the time for the diffracted intensity fall-off in Fig. 5. Possible 
inhomogeneous laser intensity distribution of both the pump and the probe beams over the 
surface induce further complication into the interpretation of the diffracted beam intensity 
behaviour. 
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Fig. 6. Schematic representation of exponential temperature decrease (top graph) with depth in 
the skin layer, and the atomic displacement A along the x-axis (lower graph) due to 
thermal pressure (negative shift toward the laser beam) and the gradient pressure (positive 
along the laser beam); l s = 28 nm for 800 nm wavelength in bismuth. 



D. The diffracted intensity decrease in the ultra-fast experiments and its relation to the classic 

Debye- Waller factor 

Many authors [12,14,15] use the classic Debye-Waller factor for determining the lattice 
temperature in a solid swiftly excited by the ultra-fast laser. Let us discuss the validity of 
classical (equilibrium) approach at the temperature approaching the melting point or in the case 
when atomic displacements are of "well-directed motion ". The decrease in the beam intensity 
diffracted from a medium where the atoms are displaced at a small (u « G A ) distances u 
expresses through atomic displacement and the reciprocal lattice vector G [24]: 

— = (eip (-/G ■ ufi = 1 - i(G ■ u) - - l(G ■ u 

7 2 \ 

The angular brackets means averaging over the distribution function. The condition 
(G-uj =0 holds if u is a random displacement uncorrected with G. That applies to the thermal 
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harmonic vibrations, («)= J uf B du, distributed in accord with the Boltzmann function 
f B ccQxp^-Mco 2 ph u 2 / k B T^, here co 2 ph is the constant phonon frequency. At the temperature 

approaching the melting point the atomic vibrations lose their harmonic character. Then, the third 
order correction to the interaction potential becomes significant, AU « Mco 2 ph u 2 -gu 3 . Now the 
average displacement from the equilibrium position (the first term in RHS of Eq.(25)) becomes 

non-zero, (uj~^—d [10,28=24] while the dependence of the average phonon energy on the 

£ b 

1 2 I i\ 3 

temperature still holds the same form as for the harmonic vibrations, -Ma> ph (u J ~ -k B T (here 

the harmonic oscillations are in three dimensions [24]). Therefore at the temperature 
approaching the melting point the classic Debye-Waller expression, ignoring the first term in 
Eq.(25), becomes invalid. The first and the second order terms of the diffracted intensity decay 
are both non-zero in respect to the averaged displacement and are proportional to the lattice 
temperature. The first order term is also non-zero in the case of coherent atomic displacement 
induced by the electronic pressure gradient. For these reasons the classic Debye-Waller 
expression cannot be used to determine the lattice temperature at the elevated temperature close 
to the melting point. 



V. DISCUSSION AND CONCLUSIONS 

We would like emphasising the several important issues missed in the interpretation of the 
experiments with Bismuth swiftly excited by ultrafast laser pulses. 

First, the recovery of the temperature dependence for the electron-phonon momentum and 
energy exchange rates, for the number density of the conductivity electrons, and thermal 
diffusion at the temperatures up to 200K above the melting point allows performing the 
quantitative interpretation of the experiments with the ultra- fast excitation of Bi. 

The number of electrons in the conduction zone from the published data [7,18-20] at the 
temperature in excess the room temperature ranges from 40% to 100% of the valence electrons. 
Modification of the potential surface by the electron excitation occurs in Bi in a way similar to 
that in other metals such as aluminium [16] or gallium [32], and can be related to the increase in 
the inter-atomic spacing due to directed, coherent displacement [9-11]. It is worth noting that 
coherent displacement is rather similar to thermal expansion. 
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Secondly, the temperature dependence of the electron-phonon energy exchange time allows 
unequivocally establish purely thermal nature of the ultra-fast laser induced transformations in 
Bismuth: the electron-phonon temperature setting up and energy equilibration occurs faster than 
the transformation takes place, that is in accord with the observed transformations in aluminium 
[16] and gallium [32]. 

Thirdly, the laser-induced inhomogeneity in the absorbed energy density and temperature 
results in the significant atomic displacements, comparable to the inter-atomic distance, might be 
a major source for the violation of the Braggs condition and the diffraction intensity decrease 
when the electron and x-ray probes were used. 

The additional concern in interpretation of laser-induced melting in Bi relates to the fact that 
liquid bismuth in equilibrium conditions is denser than solid. In the typical experimental set up 
for the ultra-fast excitation the temperature is a maximum at free vacuum-sample boundary, 
therefore the importance of expansion is obvious. A special diagnostics is needed to distinguish 
the effect of expansion on the diffracted probe from the phase transition changes in a sample. 

We also demonstrate that the use the classical Debye -Waller factor for the interpretation of 
the experiments in non-equilibrium conditions should be revised: the first order term in the 
atomic displacement at a temperature close to the melting point (that is zero for harmonic 
vibrations) is of the same order of magnitude as the second-order term. The use of the classical 
factor for the defining the temperature is questionable. 

One should note that the decrease in the phonon frequency when the temperature approaches 
the melting point signifies the onset of the vibration instability [25,26,33,34]. This is the first in 
a succession of instabilities preceding melting [34]. The frequency decrease relates to the lattice 
reconstruction and it does not indicate the termination of atomic vibrations. 

To conclude we have shown that the fast fs-laser induced excitation of Bi in all experiments 
published so far indicate the entirely thermal nature of the observed transformations. 
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